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We describe linear and nonlinear transport across a single impurity Anderson model quantum 
dot with intermediate coupling to the leads. The coupling is large enough that sequential tunneling 
processes alone do not suffice to properly describe the transport characteristics. Upon applying a 
density matrix approach, the current is expressed in terms of rates obtained by considering a very 
small class of diagrams which dress the sequential tunneling processes by charge fluctuations. We 
find zero bias anomalies of the differential conductance versus applied bias, an enhancement of the 
conductance with decreasing temperature as well as the onset of universality of the shape of the 
conductance as function of the temperature. We can address the case of a spin-degenerate level 
split energetically by a magnetic field and show that, if we assume in addition different capacitive 
couplings of the two spin-levels to the leads, one of the resonance peaks is vanishing. Furthermore, 
we get qualitative agreement with experimental data showing a cross-over from the Kondo to the 
empty orbital regime. 

PACS numbers: 73.63.-b, 73.63.Kv 



I. INTRODUCTION 

The single impurity Anderson model (SIAM) 1 has be- 
come a useful tool to describe phenomena arising in 
quantum dot devices at low temperatures. It encom- 
passes single-electron tunneling phenomena 2 , cotunnel- 
ing and resonant tunneling 3 as well as Kondo 4 physics. 
These phenomena have been verified in many experimen- 
tal quantum dot set-ups realized at the interface of a 
two-dimensional electron gas 5 10 , carbon nanotubes 11 15 
and quantum wires 16 ' 17 as well as in single-molecule 
junctions 18 . At thermal energies larger than the cou- 
pling T to the leads transport can be suppressed at low 
bias due to Coulomb blockade. In this regime sequen- 
tial tunneling dominates the transport across the SIAM, 
and tunneling in and out of the dot is well described in 
terms of rate equations 19 ' 20 , with rates obtained within a 
second order perturbation theory in the tunneling Hamil- 
tonian (i.e., first order in Y). When the temperature is 
decreased to values of the order of the tunneling cou- 
pling T or lower, the sequential tunneling approximation 
breaks down, as processes of higher order in T start to 
become important. One consequence of this is the oc- 
currence of a zero bias maximum 7 ' 21 25 or minimum 7 ' 26 
of the nonlinear conductance for small temperatures in a 
quantum dot with large Coulomb interaction, depending 
on whether the single particle resonance lies deep be- 
low (Kondo regime) or above (empty-orbital regime) the 
Fermi level, respectively. 

The intermediate coupling regime with tunnel cou- 
plings of the order of the thermal energy or higher has 
not been much investigated so far. This is in part due to 
the difficulty of developing theories capable to cope with 
strong Coulomb interactions and intermediate coupling 
at the same time. However, it is this intermediate regime 
which might be of interest for transport through some 
single- molecule junctions 27 ' 28 and is relevant to interpret 
experiments on negative tunneling magnetoresistance 12 . 



In the experiments 27 ' 28 a conductance gap is observed 
at low bias which suggests the presence of charging ef- 
fects. The gap is followed by conductance peaks whose 
broadening is larger than the estimated temperature, be- 
ing a hint that tunneling processes of high order might 
be responsible for the broadening. In Ref. 12 Coulomb os- 
cillations of the conductance versus the gate voltage are 
clearly seen; however, the occurrence of a negative mag- 
netoresistance requires the presence of level shifts due to 
higher order charge fluctuation processes 29 ' 30 . 

In this paper we describe the transport by using a di- 
agrammatic approach to the stationary reduced density 
matrix of the quantum dot and the stationary electron 
current onto one of the leads. Along the same lines as in 
Ref. 29 we include all possible diagrams which dress the 
second order tunneling rates by charge fluctuations in 
and out of the quantum dot. Different from the method 
in 29 , we do not only extract tunneling induced level shifts 
from the analytical expressions. We calculate transition 
rates and express the density matrix and the current in 
terms of those. Our "dressed second order" (DSO) dia- 
grams are a small subset of the diagrams kept within the 
so called resonant tunneling approximation (RTA), first 
proposed by 31 to describe the beyond sequential tunnel- 
ing regime. In particular, for a spinless quantum dot the 
RTA is exact and reproduces e.g. the expected Breit- 
Wigner resonance shape of the linear conductance. The 
much smaller DSO subset, too, yields the known exact 
results for the current. 

The DSO yields the rates in a straightforward way. 
We compared the predictions of the RTA and DSO both 
for the linear and nonlinear conductance in the case of 
infinitely large Coulomb-repulsion and found only small 
deviations in the intermediate coupling regime. Larger 
deviations are seen at lower temperatures where the con- 
ductance obtained by the DSO is remaining considerably 
below the RTA-result. 

In this regime both approximations yield universality 
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of the linear conductance: There is a temperature Tk 
such that the conductance is a universal function of the 
ratio T/Tk- We compared our expression for Tk with 
results of other theories for the Kondo temperature in 
the case of infinite Coulomb interaction. Interestingly, as 
in the RTA 32 , we find an exponent differing precisely by 
a factor of two from the result in Ref. 33 . 

We also address the case that the two spin levels are 
split energetically by a magnetic field and reproduce the 
result that the zero bias resonance of the conductance 
versus the bias splits up into two peaks 26 . Moreover, 
we consider the situation that in addition to an energy 
difference the levels have different capacitive couplings 
to the contacts (still equal tunneling coupling). In this 
case we obtain that one part of the double peak vanishes. 
Measurements of this kind were reported in Ref. 7 . 

Finally, we focus on the case of finite but still large 
Coulomb interaction and consider the linear conductance 
as a function of the gate voltage and the temperature. 
The effect of changing the gate voltage is a shift of the 
relative position of the level energy with respect to the 
Fermi level. This offers the possibility to investigate the 
cross-over from the Kondo regime to the mixed valence 
and finally the empty orbital regime, corresponding to 
single particle energies lying deep below, in the vicinity 
(within an uncertainty of the order of T) and above the 
Fermi level of the leads, respectively. We compare with 
experimental results in Ref. 6 and obtain in many respects 
qualitative agreement. 

The structure of the paper is as follows: Section II in- 
troduces the model of the transport current. Section III 
illustrates the diagrammatic approach and recalls known 
results for the reduced density matrix and the current 
in second order in the tunneling Hamiltonian. Analyti- 
cal expressions for the current and the reduced density 
matrix are provided in terms of rates. 

The DSO approximation is explained in Sec. IV. In 
Sec. V, VI, and VII the DSO is applied to the spinless 
case, to the case of infinite interaction and finite inter- 
action, respectively. In particular, in Sec. V and VI the 
DSO and RTA predictions are compared; the case of en- 
ergetically split levels is considered. In Sec. VII on the 
other hand we compare with the experimental results in 6 . 
Finally, conclusions are drawn in Sec. VIII. 

II. BASIC MODEL 

A. Hamiltonian 

The Hamilton operator of our system is H = Hr + 
Hq + Ht- In the reservoirs we assume noninteracting 
electrons. Correspondingly, we choose 

In this formula, the indices Z, a and k denote the lead, the 
spin and the wave vector of an electron level in the con- 



tacts, respectively; is the band energy corresponding 
to this electron level; c^k is the annihilation operator of 
the level lak and the dagger denotes the Hermitian con- 
jugate. 

The Hamiltonian of the isolated quantum dot is 

a 

where U is the Coulomb interaction and d a and S a are the 
annihilation and creation operator of the level a =t / I 
on the dot. Alternatively, the Hamiltonian of the isolated 
dot can be written as 

H Q = E o \0)(0\ +Y,E a \a){a\ + E 2 \2)(2\. 

a 

For any of the four many particle states a = 0, t, -Jo 2 
we use "E a " to denote the energy of this state. The 
Hamiltonian of the quantum dot is diagonal in the basis 
given by these four states. By comparison with the above 
representation we get: Eq = 0,£?2 = U + ^2 a E a . Only 
differences between energies of quantum dot states are 
relevant. For this we use the terminology 

E a b '•= E a — E\>. 
The tunneling Hamiltonian, 

Ht = Tike d^ciko- + H. c. (Hermitian conjugate), 

Jerk 

connects electron levels on the leads with the level on the 
quantum dot 34 . With the term "order" of a tunneling 
process/ diagram we mean its order in the "tunneling 
coupling" or "in Ht"- 

B. Initial condition 

We assume that there is an initial time at which the 
systems are still separate and express this by writing the 
initial density matrix as product of density matrices of 
the quantum dot and the leads: 

P(to) = Po(^o) ® Pr, 

where Po(to) is some arbitrary initial density matrix de- 
scribing the state of the dot; p R = p R ,i e ft ® pR.right is 
the density matrix of the leads in thermal equilibrium. 
Specifically, we choose 

where \i\ is the chemical potential of lead I and where n\ is 
a normalization factor. After this initial time we assume 
that the time evolution of the density matrix is ruled 
by the total Hamiltonian H according to the Liouville- 
von Neumann equation 35 which is the analogon of the 
Schrodinger equation for density matrices. 
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C. Thermodynamic limit and the current 

For each of the leads we define an electron counting 
operator as Ni = Xlkcr c lka c ikcr anc ^ ^he °P era tor of the 
particle current onto that lead as I\ = TV/]. Then 
the current onto the chosen lead at time t is 

jtrace (N lP (t)) = Tr (I lP (t)) =:< // > (t). 

We define the stationary current by letting the time go 
to infinity and taking the average current: 

rOO 

< Ii > oc = lim A / dt<h> (t)e- x ^- t0 \ (1) 

J to 

where A > is the argument of the Laplace transform 
of the function < I\ > (t). The total weight of the mul- 
tiplicant Ae _A ( £_£ °) over [to,oo[ is always unity, but for 
smaller and smaller values of A it will be distributed over 
a larger and larger time interval. The current in this defi- 
nition is zero as long as the contacts are finite. Therefore, 
we first let the size of the contacts go to infinity and re- 
define 

<I l >(t)= lim <J z >(t,y). (2) 

Then the current in the definition of Eq. (1) is our model 
of the dc-current measured in transport experiments with 
quantum dots. 



III. DIAGRAMMATIC APPROACH 

A. Basic method 

An analysis of the time evolution of the current, 
Eq. (2), shows that it can be separated into subse- 
quent smallest segments, so-called irreducible tunneling 
processes 32 ' 36 . The calculation of the stationary current 
can be reduced to the calculation of the corresponding 
transform of these irreducible segments. The theory is 
exact. 

The irreducible segments of the time evolution of the 
current are called "kernels". We distinguish between the 
" density matrix kernel " K which determines the reduced 
density matrix of the quantum dot and the "current ker- 
nel" Kq which defines the relation between the reduced 
density matrix and the current. The fact that the time 
evolution of the current is completely determined by the 
kernels can be expressed in a compact way by the two 
equations: 

<Ii>(t) = Tr f dt'K c (t-t')p Q (t'), 

i [ l 
Po(t) = rbo (*),#©]+ / dsK(t-s)p Q (s), 



where we use the terminology p©(t) := Tvr {p(t)} for the 
reduced density matrix of the quantum dot. The second 
equation is also called the quantum master equation 35 . 
We take the Laplace transform of both equations in the 
limit A — >> 0. Then, the second equation allows the calcu- 
lation of the stationary reduced density matrix as far as 
K(X = 0) is calculated. Finally, the Laplace transform of 
the first equation can be used to calculate the stationary 
current. 

The calculation of the current means thus the calcula- 
tion of the kernels. The contributions to the kernels are 
visualized by diagrams, their number and variety is huge. 
This forces us to take into account only special classes of 
diagrams about which we have reason to believe that they 
might be important and which we are able to calculate. 
Only in the special case of the spinless quantum dot an 
exact solution was presented 31 . 



B. Second order approximation 

Our approximation is an extension of the second order 
theory, so we explain for a moment what we mean by this 
and note its result. 

If we multiply the tunneling Hamiltonian by a 
dimension- less parameter 'V, then everything becomes 
a function of w, including the kernels u K(w) u , LL Kc(w) LL 
and finally also the current, "< // >oo (w)" . All of the 
contributions to the kernels have an order in the sense 
that the coefficients of the tunneling Hamiltonian appear 
a certain number of times. All odd orders vanish. Thus, 
they have the structure: 

K(w) = w 2 K^ + w 4 K^ + . . . , 
K c (w) = w 2 K^ +w 4 K^ + .... 

For small values of w, i. e., for weak tunneling coupling, 
one takes into account only the approximations for the 
kernels of order 2n and calculates the current on this 
basis. One obtains the current of order 2n, u I^ 2n \w) c \ 
i.e., the Taylor expansion of the current around w = of 
order 2n. The current, Eq. (1), is analytic in w 2 . In the 
case n = 1 one obtains the second order current, which 
is just of the form w 2 • constant. 



1. Second order density matrix 

One of the diagrams visualizing the second order con- 
tributions to the density matrix kernel is shown in figure 
1. A possible way of describing the process is to say that 
the quantum dot is at first in the un-occupied state "0" ; 
then an electron with spin a tunnels in two steps onto the 
dot. Finally, the dot is in the state "a". The analytical 
expressions which correspond to the diagrams are given 
by diagrammatic rules, e.g. 26 ' 31 ' 36,37 . The expression for 
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FIG. 1: An example of a second order diagram; by convention 
we let the time increase from the right to the left. The hor- 
izontal lines are called "contours", the third line "tunneling 
line" ; it represents an electron from lead / with spin "a" and 
wave vector "k" which tunnels in two steps onto the dot. The 
intersection points of the tunneling lines with the contours are 
called "vertices". The vertices separate the contours into in- 
tervals, and to each of these we assign a quantum dot state. 
The particle number of neighbouring quantum dot states can 
differ only by ±1. The chronological order of the density ma- 
trices of the quantum dot in this diagram is found as follows: 
We imagine a vertical line which cuts each of the contours 
one time. We consider then especially the intervals between 
two neighbouring vertices which are cut by the vertical line 
and take the two quantum dot states assigned to them. If 
they are "a" and "6" on the lower and on the upper contour, 
respectively, then the current quantum dot matrix is given by 
\a){b\. The imaginary vertical line we move from the right to 
the left end of the diagram. The sequence of quantum dot 
matrices in the above diagram is: |0)(0|, |0)(cr| and finally 
\a)(a\. 



the diagram in figure 1 reads: 

h h\ + i(sik a - E 10 ) ' 

where we let A, the argument of the Laplace transform, 
still be finite. For simplicity we assume degeneracy, E a = 
and write E a Q =: Eiq. Later we will consider also the 
case of different energies E a ^ E&. We let fi(e) be the 
Fermi function at chemical potential \i\ and temperature 
T, i.e., fi(s) = f((e-^)/k B T) with f(x) = 1/(1 + 
e x ). We perform then the sum with respect to the leads 
and the wave vector. The thermodynamic limit is taken 
by replacing the sum with respect to the allowed wave 
vectors k by an integral over the first Brillouin zone. The 
expression turns into: 



dkZi 



HX + z(e Zko . - E 



10J 



where Z\ is the number of allowed wave vectors in the 
first Brillouin zone per volume in the wave vector space. 
We split now the integration into two parts: First, we 
fix the band energy and integrate over the surface in the 
first Brillouin zone where £i\^ a equals this band energy. 
In a second step, we integrate over the band energies 38 . 
The integral turns into: 



de 



Me) 



hX + i(e — Eiq 



dS\T lka \' 
-e} |Ve Jff (k)|- 



There are two diagrams of second order which are con- 
tributions to the "kernel element" (a\ {K(X)\0)(0\} \a). 
They are given by the above diagram and by the one we 
get by mirroring this with respect to a horizontal axis. 
Their contributions are complex conjugate, so we have to 
take two times the real part of the above expression. In 
the limit A — » we obtain: 

(a\ {K(X = 0)|0)(0|} \a) = \ 5>+(£ 10 ), 



where we used the notation a^(e) = ai(e)fi(e) and 



ai(e) 



L 



dSZ t \T lka [ 
= e} |Ve I(T (k)| 



(3) 



In the case that the tunneling coefficients were inde- 
pendent of the wave vector the function ai(e) would just 
be proportional to the density of electron levels in lead I. 
However, we point out that we do not use such a simpli- 
fying assumption about the tunneling coefficients at this 
stage. The dimension of ai is "energy", correspondingly 
the dimension of the kernel elements is "rate" . 

The other second order kernel elements are calculated 
essentially in the same way. With the further notation 



'(e) = ai(e)(l - fi(e)) we obtain: 



{K(X = 0)\*))*\}\0) 



(2\{K(X = 0)\a)(a\}\2) 



(a\{K(X = 0)\2)(2\}\a) 



2tt 



2tt 



2tt 



i 

I 

5>r(£2i). 



For simplicity we assume for this a symmetry in the 
leads with respect to the spin, i.e., that the definition 
of ai, Eq. (3), does not depend on the spin. By con- 
sidering the contributions of diagrams as in figure 2 
one can verify that the general property of the density 
matrix kernel Tr {K\a)(a\} = holds true also within 
the second order theory. Therefore, we already calcu- 
lated implicitly the remaining kernel elements of the form 
(a\{K(X)\a)(a\}\a). 




FIG. 2: Second order diagrams of this form ensure that the 
trace of K(X)\a)(a\ is always zero. Their contributions to the 
density matrix kernel amount to — ^2 b -^ a {b\ {K(X)\a)(a\} \b). 
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For the Anderson model quantum dot, the density ma- 
trix kernel always transforms a diagonal matrix into a di- 
agonal matrix, i.e., (b\ {K(X = 0)\a)(a\} \b') = if b ^ V . 
Thus, we can say that it is a linear operator with rank 
three or lower acting on the four-dimensional space of 
the diagonal matrices since one degree of freedom is de- 
stroyed by the condition that the trace of the resulting 
matrix is zero. We can conclude that there is a diagonal 
solution "p" of the quantum master equation (QME) in 
the stationary limit. With the notations 



Paa := (a\p\a), 



1,01 



■ 1,12 



and 



2tt 



the QME turns into the following set of two equations 
(for the three variables poo,P22 and p^ = pn): 

= -poo T toi + P°° T T,oi > 
i i 

I I 

This information is sufficient to determine the stationary 
reduced density matrix since we know that the normaliza- 
tion condition, poo + 2p acr + p 22 = 1, holds. The solution 
is: 




1 Ol 1 12 
1 Ol 1 12 

V r£,r+ / 



(4) 



where we used the notations 



r ± - 

1 ab 



r ab := r+ +r 



ab * ab' 



Second order current kernel 



In order to calculate the current we have to determine 
the second order current kernel. The structure of the con- 
tributions to it is the same as that of the contributions to 
the density matrix kernel. We take into account only the 
diagrams with the final vertex on the lower contour. The 
lead-index attached to the corresponding tunneling line 
is fixed and given by the lead onto which we are calculat- 
ing the current. An additional sign, as compared to the 
density matrix kernel, has then to be taken into account. 
There are several equivalent possibilities of defining the 
current kernel 32 . 

For example, the diagram in figure 1 yields the contri- 
bution 



if 



de 



-(e) 



to the trace Tr {K c (\)\0)(0\}. 

The other contribution to this trace comes from the 
diagram in figure 2 if we set a = 0, b = a. The two 
contributions are complex conjugate and so we get: 

rr{tf c (A = 0)|0)<0|} = -2r+ 01 . 

In an analogous way one obtains: 

Tr {K C (X = 0)|2)(2|} = 2r~ 12 and 

Tr{K c (X = 0)\a)(a\} = r~ 01 - r+ 12 . 

The second order particle current is then found by ap- 
plying the current kernel to the reduced density matrix 
and taking the trace 26 : 



r toi r ^ 01 _ r ^,oi r toi 

r r,i2 r /7i2 - r /" 12 r ^i2 



,12 



(5) 

where we used "Z" to denote the opposite lead of lead "Z" 
and the abbreviations: 

TV := i7 2 r i + r+r 12 , 



Lab 



l,ab' 



The letter w denotes the coupling parameter as intro- 
duced above; for simplicity we will leave it away in the 
following. This is the particle current onto lead I. The 
net current, i.e., the sum of the two currents onto lead I 
and T, is zero. To determine the electric current one has 
to multiply by the electron charge. 

In the case of proportional tunneling coupling, i.e., 
ai = k\ol, with k\ positive scalar factors fulfilling k\ = 
1, the expression for the current can be simplified: 

{^ T I,oi ~ K i T toi) 



r(2) 



1 + Lqi _ 

r oiT 12 



+ ; 7 r 01 rr a ( KtTi > 12 ~ KlT m) ■ 



(6) 



h\ + i(e-E 10 ) 



The prefactor of the second line turns out to be the sta- 
tionary electron number on the quantum dot, i.e., the 
expectation value Tr {N^p} with 7V the particle count- 
ing operator on the quantum dot; the prefactor of the 
first line we might call the "hole number", i.e., the ex- 
pectation value of 2 — 7V . 

The second order approximation can be interpreted 
in terms of transitions. For every pair of two quantum 
dot states "a" and "6" with neighbouring particle num- 
bers and every lead I we determine a rate of transitions 
"a — » 6" caused by the tunneling of an electron from lead 
I onto the dot or from the dot onto lead Z, provided that 
the dot is in the state \a)(a\. The stationary density ma- 
trix is determined by the demand that all of these transi- 
tions compensate each other. Then we calculate the cur- 
rent by balancing the transitions. Our non-pert urbative 
approximation is an extension of the second order theory. 
The equations for the density matrix and the current in 
terms of the transition rates T^ ab still hold true but the 
expressions for these rates change. 
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IV. DRESSED SECOND ORDER DIAGRAMS 

In this section we account for diagrams similar to the 
ones of the second order theory but "dress" them by 
charge fluctuations. Figure 3 shows two possibilities of 
dressing the second order diagram of figure 1. Apart from 
the electron that tunnels in two steps onto the dot there 
is one more electron level of the leads involved. We might 
say that an electron (lower example diagram) or a hole 
(upper example diagram) tunnels for some time halfway 
onto the dot and then leaves it again. For linguistic sim- 
plicity we restrict ourselves here to speaking about par- 
ticles tunneling onto the dot, with the consequence that 
we are using the terms "hole" as well as "electron". The 
tunneling of the one electron which finally enters the dot 
is accompanied by the tunneling of further electrons or 
holes in these diagrams. 




In the first line we recognize the contribution of the sec- 
ond order diagram. However, the integrand is multiplied 
by a factor, the second line, and this reflects the par- 
ticipation of further particles. From the upper diagram 
in figure 3 we get "a~(s') L \ from the lower one we get 
"a + (£ / )". The sum yields "a(s f ) u which appears in the 
factor. 

Because of the existence of the two spins there is still 
a third way of dressing our second order diagram with 
one further tunneling line: The bubble on the lower con- 
tour of figure 3 might as well represent an electron with 
opposite spin which accompanies the tunneling. The 
contribution of this diagram is the same as the contri- 
bution of the diagram with only one spin appearing, but 
it is important because it does not have a counterpart: 
There is no way of dressing the diagram with a bubble 
on the upper contour which represents a hole of the op- 
posite spin. Finally, there is the possibility to dress the 
diagram by a tunneling line on the upper contour which 
represents an electron of the opposite spin, as shown in 
figure 4. 




FIG. 4: The existence of the state "2" leads to a fourth pos- 
sibility of dressing. The tunneling line on the upper contour 
represents an electron of the opposite spin which tunnels onto 
the dot and leaves it again. 



FIG. 3: Two examples of dressing the diagram of figure 1 with 
further tunneling lines. In the upper diagram the temporal se- 
quence of quantum dot matrices is: 1 0) (0|, |0) (a |, |0) (0|, |0) (cr\ 
and |cr)(cr|. A hole from lead l' with wave vector k' is partic- 
ipating in the process. In the lower diagram, it is a further 
electron which is accompanying the tunneling of the electron 
from the level /kcr. 

According to the diagrammatic rules the sum of the 
contributions of the diagrams in figure 3 to the density 
matrix kernel is given by 

I f d£ 

hj h\ + i(£-E 10 ) 

zl ^) 

HX + i(e - E 10 ) J h\ + i(e - e f ) ' 
where we let A still be finite and used the notations: 

i 

a := a + + a~ — ^ a\. 

i 




FIG. 5: An example of a diagram with four "bubbles" and 
one "long" tunneling line. If we count the bubbles from the 
right to the left, then for the choice of each of them we have 
essentially four different possibilities. This example-diagram 
contains all four of these possibilities. 

We saw that there are essentially four ways of dressing 
the second order diagram with one bubble. Moreover, we 
can dress the diagrams with two or even more, in general: 
n, subsequent, non-intersecting, bubbles. An example 
is sketched in figure 5. For the choice of each of these 
bubbles we have four possibilities. It might represent an 
electron of the same or the opposite spin and thus appear 
on the lower contour or represent a hole of the same spin 
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or an electron of the opposite spin and appear on the 
upper contour. The sum of the contributions of all of 
these diagrams to the kernel element (a\ {K(X)\0)(0\} \<r) 
is 

1 f d£ q + 00 fW -1 ) n 
\J V + i(g - e 7 ) 77 + i(g + g 7 - £20) J 



the time space as the Laplace back transform of this right 
hand side 39 and in this way see that its Laplace trans- 
form exists for all values of A with strictly positive real 
part. The right hand side of the above equation, too, is 
holomorphic in A on the whole half plane {Re > 0} where 
the real part is positive. Thus, we have two holomorphic 
functions on {Re > 0} which are equal on {Re > c}. The 
theory of holomorphic functions says that they must be 
equal everywhere. 



/ 



de 



a+(g)/fi 



rj+i(e-e') r)+i(e+E' -E 2 q) 



where we replaced "ftA" by "77" . This summation of a ge- 
ometric series can be justified for "large enough" values 
of A. However, we are interested in the limit A — » 0. For 
this we remember that we want to calculate the Laplace- 
transform of the sum of the corresponding diagrams in 
the time space and let A be complex. If the real part 
of A is sufficiently large, say, ReX > c, then the Laplace 
transform of the sum of the diagrams in the time space 
is indeed given by the right hand side of the above equa- 
tion. One can represent now the sum of the diagrams in 



A. DSO tunneling rates 



All of the other second order diagrams can be dressed 
in the same way. For the diagrams connecting the par- 
ticle numbers one and two we see that a support of the 
tunneling by holes and electrons of the same spin and 
by a hole of the opposite spin is possible, but not by an 
electron of the opposite spin. We obtain the following 
transition rates within this dressed second order approx- 
imation: 



2n r & af(e) [(a + q+)(g) + a+(E 20 - e)} 

h J 8 : 7r 2 [(a + q+)(g) + a+(E 20 - e)f + [g + p a + a+ (s) - p a+ (E 20 - e) - E 10 } 2 ' 

2tt f af(e) [(a + q~)(g) + a~(E 20 - e)} 

fi J £ n^(a^a-)(e)^a-(E 20 -e)} 2 ^[e^p a+a -(e)-p a -(E 20 -e)-E 21 ] 2 ' 



where we define for any function h the function ph by 

p h (e) := / duj^ '——i (9) 

In the case of weak coupling these transition rates in- 
deed turn into the transition rates of the second order 
theory. More precisely, we get if we multiply the tunnel- 
ing Hamiltonian with a parameter w as done in IIIB: 

r^ ab (dressed second order) (w) 

— ■ x — — — > l(w 0). 

ab (second order) (w) 

The stationary reduced density matrix within the DSO 
is given by Eq. (4) and the current is given by Eq. (5) 
and in the case of proportional coupling by Eq. (6). 

In Ref. 26 a diagram selection called the "resonant tun- 
neling approximation" (RTA) was applied: The case of 
infinite interaction was considered and therefore the di- 
agrams containing the state "two" do not contribute; 
transition rates were derived, the stationary density ma- 
trix was determined and finally the current was obtained. 



The RTA takes into account all of the diagrams within 
the DSO which do not contain the state "2". However, 
there are many diagrams outside the DSO which are con- 
tained in the RTA as we will discuss in Sec. V. Compar- 
isons between the predictions of the DSO and RTA will 
be performed in the cases of the spinless quantum dot 
and of the SIAM with infinite interaction. 



B. Linear conductance within the DSO 

We assume that we can obtain the second order func- 
tions ai by multiplication of the density of electron levels 
in the leads by a coupling constant while in general their 
definition is more complicated (Eq. (3)). About the den- 
sity of electron levels we make simplifying assumptions 
such that we can concentrate on effects which are not 
due to special behaviour of the density of electron levels. 
In particular, we place the Fermi level in a point with 
respect to which the density of electron levels is symmet- 
ric. Then the chemical potential at equilibrium coincides 
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with the Fermi level for all temperatures, see figure 7. respect to the bias, eVuas = Pi — Pi > one obtains at zero 
By deriving the formula for the current, Eq. (6), with bias the following expression for the linear conductance: 



q(DSO) = Aki 



(2-n©; 



\V bias =0 



n© 



v bias =o 



f f d£ 7T 2 a(s)[(a+a+)(s)+a+ (E 20 -s)} 



r , 7r 2 a(£)[(a+a-)(£)+q-(£;20-£)] 



V bias =0 

v bias =o 



-1 f/ f \ 

k B TJ \ k B T J 

fc s T^ ^ /c b t y y 



(10) 



where ^ is the Fermi level, f(x) 
malized Fermi function, 



n© 



1/(1 + e x ) the nor- 



r+r 12 



is the particle number on the dot as noted above, and 
where we use the abbreviations "doi(^)" and "di2(s) LL 
for the denominators in the expressions for the transition 
rates, Eqs. (7) and (8), respectively. The prefactor 4k^- 
is one in the case of symmetric coupling and less than 
one otherwise. We finally included the electron charge 
into the formula. 



V. RESULT OF THE DSO FOR A SPINLESS 
QUANTUM DOT 

We consider here the case of a spinless quantum dot 
with only two possible states "0" and "a", so with only 
one spin. One obtains it by not performing the sum with 
respect to the spin in the Hamiltonian. This problem is 
equivalent to the SIAM with = E^ and U = in the 
sense that current across the SIAM quantum dot is then 
just two times the current across a spinless quantum dot. 

The contribution of all diagrams outside the RTA to 
the kernels is zero 40 for the spinless quantum dot and 
thus, the RTA is exact in this case. The relation between 
our approximation and the diagrams of the RTA is de- 
scribed by figure 6. In the case of the spinless quantum 
dot treated here as well as in the case of infinite inter- 
action the figure characterizes the relation between the 
diagram selections completely. 

The DSO approximation concentrates on the second 
order diagrams dressed by further tunneling lines, the 
RTA takes these ones in order to construct even more di- 
agrams. The resulting diagrams are a combination of an 
integer number of DSO-diagrams. The transition rates 
of the DSO in the case of the spinless quantum dot we 
get by taking into account only those diagrams within 
the DSO which contain only the quantum dot states "0" 
and V: 
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FIG. 6: The structure of the diagrams within the RTA. When 
projecting all of the tunneling lines onto one contour only, the 
first line shows the diagrams within the DSO; in the following 
lines, two, three, in general: n diagrams of this kind are com- 
bined to create new ones out of them. For the sum of all of 
these diagrams, an integral equation was derived 31 . Its origin 
is the relation between the contributions of all combinations 
of n or less diagrams of the DSO on the one hand and of 
all combinations of n + 1 or less DSO-diagrams on the other 
hand. In the limit n —> oo the two are equal. 



We get the density matrix 



Poo 
Pn 



r + ! 



■ 01 



(?01 



and the particle current in the case of proportional cou- 
pling: 



tDSO 



de- 



a\e) (f T -fi)(s) 



7r 2 a 2 (s) + (s+p a (s)-E 10 ) 



2 ' 



In the case of proportional tunneling coupling the result 
of the DSO for the spinless quantum dot is actually the 
same which was presented within the RTA 31 and thus 
exact. In the case of non-proportional tunneling coupling 
the results become different. 

The diagram selection defined by the first line in fig- 
ure 6 we might call the "simple" selection since the pair 
formation of 2n subsequent times is one of the simplest 
possible irreducible pair formations. In the case of finite 
interaction the DSO diagram selection contains less dia- 
grams than the simple selection. The DSO for finite U 
does not describe the noninteracting limit, where U = 0, 
correctly. This can be seen by comparing the formulas 
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for the linear conductances of the DSO in case U = on 
the one hand and of the DSO in the spinless case on the 
other hand. However, we know by now that the simple 
diagram selection does describe the noninteracting limit 
correctly. Moreover, the simple selection seems to sug- 
gest the natural way to extend the RTA to the case of 
finite U, even though the diagram summation might be 
technically difficult. The simple selection we want to dis- 
cuss elsewhere. We will now concentrate on applying the 
DSO approximation to cases with nonzero interaction. 

VI. THE CASE OF INFINITE INTERACTION; 
COMPARISON OF RTA AND DSO 



( e +Pa+a+( £ ) — Eio) 2 f° r the common denominator. Dif- 
ferences are only found in the prefactor and in the nu- 
merator. To compute the conductances we have to make 
a choice about the second order function a(e). We wrote 
a(e) = a(Ep) • b(e) with a dimensionless function b(e) 
fulfilling b(Ep) = 1. The variable a(Ep) is then our 
coupling parameter. Figure 7 shows how we chose the 
function b(e). In figure 8 we compare the linear conduc- 
tances as a function of the gate voltage obtained within 
the RTA and the DSO for various temperatures. Notice 
that the shape of the curve saturates at low temperatures 
both within the RTA and the DSO; the effect of decreas- 
ing the temeperature further and further is only a shift of 
the graph. As we will show in the next subsection, these 
features do not depend on the choice of b(s). 




FIG. 7: Energy dependence of the dimensionless function 
b(e) — ol(£)/cl(Ef). We placed the Fermi level Ef in the 
middle in order to ensure that the chemical potential at equi- 
librium always equals the Fermi level. Moreover, we chose 
W — leV, w = 0.9W. A cut-off is needed in order to ensure 
the existence of the principal parts p a ,p a + etc. of Eq. (9). 

In the case of infinite interaction, U = oo, it makes 
sense to neglect all of the diagrams which contain the 
state "2", assuming poo + J2 a P™ — 1- This was 
done within the RTA and a result for the current was 
derived 26 . One can do the same with the DSO. The for- 
mulas for the linear conductance of the RTA and DSO in 
the n U = oc case" read: 




E - E c (meV) 

■-10 F v ' 



FIG. 8: Linear conductance of the DSO as a function of the 
energy difference Eio for different temperatures. The dashed 
lines show the result of the RTA; the coupling we chose to be 
ol(Ef) — 0.042meV and W — leV. (The same choices we 
made for a later comparison with an experiment, discussed in 
Sec. VII B.) A temperature of IK corresponds to a thermal 
energy of about ksT w O.lmeV. For large temperatures the 
resonance is smeared out, its centre is found roughly around 
the Fermi level. For decreased temperatures the position of 
the maximum is shifted and the width of the peak is propor- 
tional to a(Ep)- While the logarithmic shift of the peak with 
the temperature does not stop, the shape of the curve and 
its maximum value saturate. We find numerically within the 
RTA the maximum « 0.99e 2 /^, within the DSO « 0.69e 2 /h. 
We emphasize that neither these maximum values nor the 
shapes of the curves in the limit of small temperatures de- 
pend on the the way in which the function b(s) is chosen. 
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DSO 



de 
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d(e) k B T J V k B T 



4«z« r — (2-ra©) 



E f 



I 



de- 



[a(a + a + )] (e) 



d(e) 



k B T 



r 



k B T 



where we used the abbreviation d(e) := n (a + a + ) (e 



A. Derivation of a formula for the Kondo 
temperature in the infinite U-case 

For a fixed value of the gate voltage, i.e., for fixed £io, 
the linear conductance becomes a function "G(T)" of the 
temperature. This is expected to display universality 41 in 
the following sense: There is a temperature Tk such that 
G(T)/G m ax is a universal function of the ratio T/Tk, 
where G max is the maximum value of the conductance. 
This statement we can show about G DSO as well as about 



10 



qRTA Both the RTA and the £)gQ however, do not 
yield the expected convergence of G(T) — > G ma x = 4«/Kp 
2e 2 /h (T^O). 

As we will show in the appendix one can represent 
G DSO in the form 

where the definition of F DSO is 

7T 2 /W + /(*)) 



F °so M= f dx ^ip 

J TT 2 (l + f(X 



)) 2 + 



and where 



E 1Q = E F + a(E F )p t 



(E F ). 



Finally, the function (j) a ,b(x) = g(x) — a + xb + log(b) is de- 
fined by the use of Eq. (15). The temperature T a ^ Ep ^ we 
define by the demand k B T a ^ Er ^ = a(E F ). The function 
b+(s) is given by 6+(e) = b(e)f((e - E F )k B T). 

The particle number n is a function of the tunnel- 
ing rates and still contains the temperature. However, 
for temperatures k B T <C a(E F ) these become essentially 
independent of the temperature such that we can con- 
centrate on the temperature dependence of the rest. The 
integral with respect to x contains the derivative of the 
Fermi function and is thus concentrated in a region of the 
order of one around zero. Therefore, we can in the case of 
small temperatures compared to a(E F ), k B T <C a(E F ), 
neglect the linear term in 4> a ,b(x) and estimate 

<t>a,b(x) &g(x) -a + log(b). 

The simplification enables us to write 



G 



DSO 



4/^-— (2 - n ) 



7 DSO 



E 



10 



E 



10 



log; 



with 



F 



DSO 



J TT 2 (l + f(x 



a{E F ) J-a(E F ) 

^f'{ x){l + f{x)) 
W + (g(x) + c) 



2 ' 



For large positive values of c, F DSO (c) takes small posi- 
tive values. At some value "c max " a maximum is reached. 
In between there is a value of c, "ci/ 2 ", where we have 

F DSO {ci/2)=05F DSO {Cmax y 

We define the Kondo temperature by the demand c = 
ci/2, i.e., 



1 K 



a C l/2 



exp 



^10 



E f 



a(E F ) 



T a (E F ) exp -p b + (E F ) 



The second line seems to depend on the coupling a(E F ), 
but this dependence is weak because of the logarithmic 
dependence of p b +(E F ) on the temperature which we 
show in the appendix. The dependence on the band- 
width W as introduced in figure 7 is proportionality as 
long as a(E F ) <C W such that we arrive at 



k B T K = 7Wexp 



Eio — E F 
a(E F ) 



(11) 



after numerical evaluation of the constants. The pref ac- 
tor, in our case "7", changes if shape of the band (fig- 
ure 7) is chosen in a different way, for example, to be 
Lorentzian. The rest of the derivation is independent of 
the choice of the coupling function. 
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FIG. 9: A comparison of the universal function of the DSO 
with an NRG-fit 6 : Both functions take the value 0.5 at 
T — Tk and are normalized in such a way that the maximum 
is one. A striking feature of the fit is that within one power of 
ten the linear conductance goes up from 50 % to about 95 % 
of its maximum while our function is growing much less in this 
interval. We included also the universal curve which one gets 
for the RTA. In agreement with figure 8 and unlike the NRG- 
fit, the universal functions of RTA and DSO are going down 
again for even smaller exponents than shown here. The equa- 
tion for the NRG-fit is G(T) = G max (1/ [(T/T' K ) 2 + l]) s 
with T' K = T K /(2 1/s - 1) 1/2 so that G(T K ) = G m a X /2; we 
chose s = 0.2. 

We summarize the conditions for which the done sim- 
plifications are valid: 



k B T < a(E F ) < W. 



Then we can represent 

^2 



G 



DSO 



4KiK T —(2-n Q )F 



DSO 



Cl/2 



los £)' 



where Tk is given by Eq. (11). The linear conductance 
becomes a universal function of T /Tk] at T = T^, G(T) 
reaches one half of its maximum. In figure 9 we compare 
our result for the universal function with a fit to the one 
obtained by NRG-calculations 41 . Essentially the same 
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arguments one can apply to the RTA in order to obtain 
an analogous universality; the formula for Tk deviates 
only in the prefactor. The relation between a(Ep) and 
the coupling parameter 'T" by the use of which Tk is 
most frequently expressed, e.g. 6 ' 24 ' 32 , is T = 2iia(EF). 



B. Zero bias anomalies of the differential 
conductance 

In addition to the linear conductance we considered 
the differential conductance obtained within the infinite- 



U DSO. We notice that in qualitatively the same way 
as the RTA 26 the approximation produces a zero bias 
maximum of the differential conductance in case E\q lies 
below the Fermi level (figure 10) and a minimum in case 
it lies above or in the vicinity of the Fermi level. The 
effect is getting more pronounced for smaller and smaller 
temperatures. 



The generalization of the infinite-U DSO to the case 
of different energies E^ ^ E^ is straightforward. One 
obtains the tunneling rates 



J 



1 la 



_ 2lT f 



de- 



*?(e) (a(e)+a+(e + E Str )) 



tt 2 (a(e) + a+{e + E Sa )y + (e - E a0 + p a (e) + p a+ (e + E 9ar )y 



(12) 




eV bias (meV) 



FIG. 10: The differential conductance versus the bias. We set 
E$ =E F - lmeV and chose a(E F ) = 0M2meV, W = leV 
as in figure 8. We fixed one of the chemical potentials at the 
Fermi level, /j,i = Ef, changed only the other one and de- 
fined eVbias = fjif — E F - We assumed a capacitive coupling 
between the leads and the quantum dot in such a way that 
Eio(Vbias) — e[°q + 0.2eVbias- We see a resonance appear- 
ing at zero bias for small temperatures. The dashed lines 
show the result of the RTA. The resonance is becoming more 
and more pronounced with decreasing temperature as it is 
observed in experiments 10 . The shape of the curve depends 
on the capacitive coupling and on how the window between 
the two chemical potentials is opened; however, the appear- 
ing of the zero bias anomaly does not in principle depend on 
these choices as one can conclude from the fact that they are 
irrelevant for the differential conductance at zero bias. 



We see (figure 11) that the zero bias anomaly is split 
according to eVbias — =t^U' m agreement with theo- 
retical results 23 ' 26 and experiments 7 ' 25 . Small deviations 
from this behaviour were observed in Ref. 17 . 

In Ref. 7 a resonance close to zero bias whose po- 
sition changed slightly with the gate voltage was re- 
ported. The dependence of the position on the gate volt- 
age was explained by the conjecture that two different 
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FIG. 11: The zero bias resonance is split up if E^ ^ E±. 
We chose the temperature T — lOOmK and the remaining 
parameters as in figure 10 apart from the splitting, E a o = 
E F - lmeV ± E a& /2. 



wave- functions (not only two different spins) might be in- 
volved, such that the assumption of different capacitive 
couplings of the levels to the gate electrode was justi- 
fied. However, in this case one would expect to see a sec- 
ond peak at the opposite bias. This was not measured. 
We assumed different capacitive couplings of the levels 
to the leads and obtained that with growing asymmetry 
one of the peaks is changing position, getting wider and 
much less pronounced. The other one, however, is get- 
ting sharper (figure 12). The explanation for this at the 
level of the transition rates, Eq. (12), is that T^ changes 
considerably with the bias in regions where 

fii — Hi = E ad - 

since then the region of high values of p a +(e — E^ a ) is 
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FIG. 12: An asymmetry of the capacitive couplings of the 
levels to the leads has the effect that the resonances become 
asymmetric; one is starting to vanish, the other one is getting 
sharper. The black line is the one already appearing in figure 
11. For the red line we changed the capacitive couplings in 
such a way that E^Vbias) = E$ + OAeVbias, E^Vbiaa) = 



E ( ^ +0.1eV bi as 



leaving or entering the interval over which the integral 
essentially goes. This leads to the condition "eVbias — 
±-Et4-" f° r ra pid change of the current with the bias. In 
case of different capacitive couplings of the levels to the 
leads the energy difference becomes a function of the bias. 
With increasing bias, one of the differences is decreasing 
while the other one is increasing. Thus, one of the reso- 
nances is getting sharper while the other one is smeared 
out. The positions are no more symmetric with respect 
to zero bias. 

Moreover, we noticed that also asymmetric tunnel cou- 
pling can have the effect that one of the resonances is 
getting less pronounced. One can let the coupling func- 
tions oli(e), Eq. (3), be dependent on the spin as well 
as on the lead and thus obtain further independent pa- 
rameters. This we did and we can qualitatively confirm 
the assumption that different tunnel couplings of the lev- 
els to source and drain, too, can be responsible for the 
observation of only one peak 7 . 



VII. LINEAR CONDUCTANCE AT FINITE U 

In this section we investigate the linear conductance 
at finite interaction according to Eq. (10). By consid- 
ering the integrals we can see qualitatively that we can 
expect an enhancement of the conductance with decreas- 
ing temperature if Eiq lies below and E 2 i lies above the 

Fermi level: The function /' ( jk^f ) °^ £ ^ as total 
weight one and is concentrated in a region of the size of 
the thermal energy around the Fermi level. In the de- 
nominators (Eq. (7)) the behaviour of the p- functions 
becomes important. The principal part Ph(z), Eq. (9), 
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13: Behaviour of the normalized function 6^ := 



a + (e)/a(EF) for which we use the abbreviation and the 
resulting behaviour of the function p b + (e) := f du)(b^ (e+cj) — 

b^(e — (jj))/uj around the Fermi level. It measures the surplus 
of electrons with energy larger than e with a strong emphasis 
on the situation locally around e. Therefore, the principal 
parts display a dip around Ef- If we decrease the tempera- 
ture by a factor of ten, then the value of the principal part in 
the centre goes down by the logarithm of ten. Moreover, if we 
stretch the narrower of the two dips by this factor, then we 
obtain the other dip. There is a well defined universal shape 
of the dips as we will show in the appendix. 



measures an asymmetry of the function h with respect 
to e (h is an arbitrary function here). In particular, p a + 
takes negative values around the Fermi level since in this 
region a+(e) is decreasing (figure 13). 

If we decrease the temperature, then the decay of the 
values of a + will be more rapid and thus the absolute val- 
ues of p a + (e) are getting larger. Around the Fermi level, 
e + p a + (e) approaches the energy difference Eiq and the 
integral increases. At some point the sum even reaches 
and crosses the level position and then the integral de- 
creases again. For the other integral the arguments are 
analogous. The energy correction p a - is positive here 
and increases if we decrease the temperature. 



A. DSO-Conductance for weak and strong coupling 

For the numerical implementation of Eq. (10) we wrote 
the coupling functions a(e) still as a(Ep)b(e)^ where the 
choice of b(e) is given in figure 7. We modeled the cases 
of strong and weak tunnel coupling by large and small 
factors a(Ep), respectively. In figure 14 we show plots 
of G DSO (E10) for different values of the tunnel coupling. 
In the limit of weak coupling we reproduce the result 
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of the second order theory while for increased coupling 
we expect essentially three effects: The peaks are getting 
higher, broader and the maxima are moving towards each 
other. 
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FIG. 14: Plots of the differential conductance at zero bias 
as a function of Eio (corresponding to a plot as a function 
of the gate voltage) for different tunnel couplings. We chose 
the interaction to be U = lmeV and the temperature T = 
lOOmK. For weak coupling we see peaks of small height whose 
width is given by the temperature; the peak positions are 
quite precisely defined by the resonance conditions Eio = Ef 
and E21 = Ef- For increased coupling the corrections are 
becoming more and more important; the width is increasing 
with C6(Ef), the peak position is shifted by the corrections 
p a ± . The DSO approximation breaks for strong couplings 
where it produces a strange sharp dip in the centre of this 
plot which is not observed in reality. 

For strong coupling and for values of Eio ~ Ep — U/2 
we expect irreducible tunneling processes outside the 
DSO-approximation to become more relevant: The pop- 
ulation of all of the four possible states of the quantum 
dot is of the order of one here. Hence, processes which 
transfer |0)(0| into |2)(2| or vice versa or |cr)(<j| into |<t)((t| 
should become more relevant. This might be an explana- 
tion for the failure of our diagram selection in describing 
this regime of parameters. 



B. Prom the empty orbital regime to the Kondo 
regime 

We want to compare now the result of the DSO for 
the linear conductance as a function of Eio and of the 
temperature with experimental data 6 . In the experiment 
a region within a two dimensional electron gas was iso- 
lated by electrostatically generated tunneling barriers. In 
this way a quantum dot which is tunnel coupled to leads 
was formed. The linear conductance was measured as a 
function of the gate voltage and the temperature. The re- 
sults were interpreted in terms of the SI AM. The authors 
distinguish between three different regimes of parame- 
ters, depending on whether the level position, i.e., Eio, 
is far below the Fermi energy (here the particle number is 
one, "Kondo regime"), in the vicinity of the Fermi level 



("mixed valence regime") or above the Fermi level (here 
the particle number is zero, "empty orbital regime"). We 
tested the performance of the DSO approximation by ad- 
justing the parameters and comparing with the experi- 
mental data. 
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FIG. 15: A plot of the full widths at half maximum 
("FWHM") of the peaks in figure 16 as a function of the 
temperature; the value of a(Ep) is here 0.042meV . For 
small temperatures the FWHM seems to saturate at a value 
of about 0.39meV which is in agreement with the value of the 
experiment 6 . For large temperatures the FWHM increases; 
the graph has positive curvature which is happening since the 
two peaks are getting mixed as we increase the temperature. 
This is why we cannot define the coupling constant in the 
same way as done in the experiment: There a linear depen- 
dence of the FWHM on the temperature for large tempera- 
tures was observed and the coupling constant was determined 
via the slope of the plot. 

We fitted the parameters in the following way: The 
temperatures are given explicitly. The value of the 
Coulomb interaction U, too, we take directly from the 
experiment. To fix the coupling parameter a(E F ) we 
plotted G(Eio) for various values of it. We determined 
the parameter by the demand that the full width at half 
maximum of the peaks is close to the measured values 
(figure 15). 

In figure 16 we show a plot of the linear conductance as 
a function of the energy difference Eio for different tem- 
peratures. We get qualitatively very similar behaviour 
as in 6 . With decreasing temperature, the peaks are mov- 
ing towards each other, they are getting higher and their 
widths are getting smaller and seem to saturate finally. In 
the end, we adjusted also the factor Akikj which expresses 
an asymmetry of the tunneling couplings to source and 
drain. This, however, is only a fit made in such a way that 
the absolute values of the linear conductance are about 
the same in theory and experiment. The experiment was 
addressed also by Ref. 42 where a different asymmetry was 
assumed and good quantitative agreement was obtained. 
Therefore, it is difficult to formulate exact rules defining 
in what regime of parameters the DSO is quantitatively 
correct. 

For a further comparison we show the dependence of G 
on the temperature for fixed Eio- In order to do this we 
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FIG. 16: A plot of the linear conductance as a function of the 
energy difference Eio which corresponds to a plot as a function 
of the gate voltage. We chose the interaction U = 1.9meV, 
the coupling q>(Ef) = 0.042raeV, and the temperatures in 
agreement with the experiment 6 ; the asymmetry of the tunnel 
couplings to left and right lead we chose as Akikj = 0.5, i.e., 
we assume an asymmetry of about ki : Kf = 0.17. For even 
smaller temperatures we would again see the sharp dip in the 
centre of the plots already shown in figure 14. 



express E\q in terms of its position relative to the Fermi 
level and divide it by some quantity T which characterizes 
the tunnel coupling in a similar way as in the experiment. 
By the plot we conclude that we have agreement of our 
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log(T/mK)/log(10) 

FIG. 17: A plot of G(T) for fixed values of the gate volt- 
age, expressed in terms of e~o := (Eio — Ef)/T. Here we take 
r « 0.3meV as in the experiment with the argument that the 
saturation width of figure 15 is the same as in the experiment. 
The different parameter regimes are called "empty orbital 
regime" (eb > 0), "mixed valence regime" (—0.5 < Sq < 0) 
and "Kondo regime" (so <C —0.5). 

theoretical result with the experimental data in the sense 
that the transition from the empty orbital regime, where 
with decreasing temperature we see only a decrease and 
then rather constant behaviour, to the Kondo-regime, 
where we see only an increase, happens within an interval 
of £?io-values of the size of about 2T. 

On the other hand, we see a qualitative deviation in 



the regime Eiq « Ep: G DSO displays quite clearly a 
decrease with decreasing temperatures for small tem- 
peratures while in the experiment this decrease is much 
weaker. 



VIII. CONCLUSIONS 

We used a diagrammatic approach in order to describe 
transport across a SI AM quantum dot. We found a 
minimum selection of diagrams which we called " dressed 
second order" (DSO) diagrams, which straightforwardly 
yield the current in terms of transition rates. In general, 
the DSO represents the natural extension of the sequen- 
tial tunneling approximation, valid for large interactions 
and when ksT ^> T, to the regime ksT ~ T. In partic- 
ular, the conductance versus gate voltage exhibits peaks 
with a broadening no more given by the temperature but 
by the tunnel coupling. 

Moreover, the diagram selection contains a zero bias 
anomaly developing at low temperatures. We showed 
that, if the degenerate level lies below the Fermi energy, 
then it is a zero bias maximum of the differential differ- 
ential conductance which appears for low temperatures 
and is getting more pronounced if the temperature is de- 
creased further. This is in qualitatively good agreement 
with experiments 7 ' 10 . 

We showed that the anomaly displays features of the 
Kondo effect such as a universality in the dependence of 
the linear conductance as function of the temperature. 
We investigated the behaviour of the anomaly in case 
a magnetic field is applied and discussed the impact of 
asymmetries with respect to capacitive or tunnel cou- 
plings to the leads. 

We pointed out the close relation of the DSO diagram 
selection to that of the resonant tunneling approximation 
(RTA) and compared their results in the case of infinitely 
large interaction. The RTA is more precise in the sense 
that it includes more diagrams; however, the DSO can 
be applied more easily to the case of energetically split 
levels. 

Finally, we applied the approximation in the case of 
finite interaction and compared its result at the level of 
the linear conductance with an experiment, Ref. 6 . We 
found good qualitative agreement with the experimen- 
tal data. It is difficult to formulate general rules defining 
for what regime of parameters the DSO approximation is 
quantitatively correct. Being an extension of the sequen- 
tial tunneling approximation, we expect the DSO to well 
cover the regime of large interactions down to temper- 
atures fc#T ~ T. Obvious failure of the approximation 
in describing the experiment 6 we observed in the region 
Eio ~ l/2(Ep + Ep — U) for small temperatures. 

The DSO for finite U does not produce the observed 
plateau of the linear conductance as a function of the gate 
voltage for small temperatures forming between the two 
resonances 10 . Additionally, the DSO for finite U does 
not correctly describe the noninteracting limit. There is, 
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however, a natural extension of the DSO which is indeed 
doing so, as outlined in Sec. V. Thus, there is hope that 
the DSO can be improved in such a way that the case of 
small interaction is described better. Moreover, the same 
class of diagrams could provide a natural way to extend 
the RTA to finite interaction as described by figure 6. 
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Appendix 



We want to derive a representation 



/ 



de 



d{e) 



v bias =o 



? RTA 



E 



k B T 

k B T 

a(E F ) ' a(E F ) 



1 f , 



E F 



k B T 



10 



E 



10 



where the denominator of the integral is given by d(e) := 
7r 2 (a + a + ) 2 (e) + (e + p a +a+_( £ ) ~ ^io) 2 , the function 
jpRta ^ g un i versa ] anc [ w here E\q does not depend on the 
gate voltage or on the temperature. To this end we write 
the second order function a(e) as a(E F )b(e) and divide 
numerator and denominator of the integral by a(E F ). 
Moreover, we write e = E F + xk B T and integrate with 
respect to x instead of e. We argue then that the integral 
is concentrated in a region of a few multiples of the ther- 
mal energy around E F and that it is, because of this, 
for sufficiently small temperatures allowed to estimate 
b(e) = 1 and pb(s) = Pb(E F ); pb(E F ) is zero because we 
chose the function b to be symmetric around E F . After 
these modifications we obtain the integral: 



/ 



dx- 



ir 2 f'(x) 



'7T 2 (1 + f{x)) 2 + 4> 2 { X y 

where we used the abbreviation 
k B T Ep — E\q 



(13) 



(f)(x) := x 



a(E F ) a(E F ) 



p b +(E F + xk B T) 



and where the function b J is given by bj,(e) = 
b(e)f((e-E F )/k B T). 

Thus, what remains to be done is the analysis of 
p b +(E F + xk B T). First of all, we consider its values in 
x = for different temperatures. By taking the deriva- 
tive with respect to the temperature one obtains: 

d_ KEf + *)/ (gr) ~ KEf ~ u)f (g) 
dT J uj 

i r°° 

= 7F dx ~ /'(*) KEf + xk B T) + b(E F - xk B T)} 
1 Jo 
1 

T' 



Jo 



y) 



f (b(E F + k B T(x + y))f(x- 
—b{Ep + k B Ty)f{y) 
-b(E F + k B T(x - y))f(x - y) 
+b(E F - k B Ty)f(-y)). 



For every single value of the limit T — » can be taken. 
We can guess that the limit of the integral is given by 
the integral of the point- wise limit, 



9(x) := f° 



dy 



-fix 



y) - f(y) 
-y) + f(-y))- 



(15) 



Then we replace the function of x given by Eq. (14) by 
g(x) with the argument that for small temperatures the 
deviations between the two can be expected to be small. 

For the proof we want to apply Lebesgue's convergence 
theorem, so we need an integrable upper bound which 
is independent of the temperature. Moreover, we write 

Jo°° = Jo + Ii° smce the different intervals make differ- 
ent treatment necessary. 
The integrand has the form 

(AB)(x + y)- (AB)(y) - (AB)(x - y) + (AB)(-y). 

As to the integral J^, we group the terms with equal 
"x" into pairs and consider the two resulting differences 
separately. By adding and subtracting the mixed terms 
A(x + y)B(x — y) one can see that we have even a con- 
stant upper bound within this interval. The conditions 
which we demand from the function b(e) for this are the 
following: 

• It is bounded, \b(e)\ < B, B independent of e. 

• It satisfies a Lipschitz condition of the form: 

| b(e) — b(e')\ < L\e — e'\, L independent of e and 



We mention that the Fermi function, too, has the two 
properties; the latter can be seen by using the fact that 
the derivative of the Fermi function is bounded and the 
mean value theorem. Moreover, a Lorentzian or our 
choice of the function b (figure 7) fulfills these conditions. 

As to the integral J-j 00 , we group the terms with equal 
sign in front of y into pairs. Again, it is useful to add 
and subtract the mixed terms, e. g., A(x + y)B(y). We 
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obtain then: 
1 



y 



b(E F + k B T(x + y)) (f(x + y) - f(y)) + 



- (b(E F + k B T(x + ^/)) - 6(£ F + k B Ty)) f(y). (16) 

The point-wise limit of the first line is 
1/y (f( x + y) ~ f(y))i an d the convergence is bounded 
by B /y \f(x + y) — f(y)\. This is integrable because we 
can estimate: 



\f(x + y)-f(y)\ 



< max {\f (z) | : \z 
=: m x (y) \x\ , 



■y\ < \x\}\x\ 



where we used y > 1. We treat x as a constant during 
these considerations. Because of the rapid decay of the 
derivative of the Fermi function m x (y) is integrable. 

The point-wise limit of the second line is zero. In or- 
der to get an integrable upper bound we introduce the 
function 



L>(e) 



sup < 



W)-KOI 



e" - e' 



: e <e' <e 



•}• 



Then we can estimate the second line by 
1 



< 



(b(E F + k B T(x + 2/)) - 6(£ F + k B Ty)) f(y) 
[E F + /c B T(?/ -\x\)-E F + k B T\x\] 



L>(E F 



k B T(y 



where we multiplied and divided everything by k B T \x\ y. 
(The square bracket is just a complicated way of writing 
LL k B Ty LL .) We note now that for sure L > (s) < L and 
make the further assumption that the function £L > (e) is 
bounded over any interval which has a lower bound, i.e., 
for any e$ we have 

L^o) := sixp{|eL>(e)| : £ > e } < oo. (17) 

This assumption is fulfilled both for Lorentzian shapes of 
b(e) and for our way of choosing the second order func- 
tion (figure 7), the reason being the rapid decay of the 
derivatives of these functions. Using these properties we 
obtain for temperatures smaller than some arbitrary tem- 
perature Tq the upper bound: 



We showed the convergence 

p b + (E F + k B Tx) - p b + (E F ) -+ g(x) (T -+ 0) 
for arbitrary x, where the limit is given by the definition 




FIG. 18: A plot of the function g(x) defined by eq. (15), 
i.e., the normalized shape of the functions p h +(e) around the 
Fermi level in units of the thermal energy. The growth of 
g(x) is logarithmic in the sense that xg \x) — > 1 (\x\ — > 
oo). However, because of the presence of the derivative of the 
Fermi function in the integrals, the behaviour of g(x) around 
x — is more important for us. 

(15). In figure 18 we plotted the function g(x). Already 
earlier we noticed that the dependence of p b + (E F ) on the 
temperature is logarithmic. By putting the two pieces of 
information together we can estimate: 



p b+ (E F + xk B T) w p b+ (E F ) + log 



L a(E F 



+ g{x), 



where T a ^ Ep ^ is defined by the condition k B T a ^ Ep ^ = 
a(E F ). We insert this into the integral (13) and obtain: 



(p(x) 



T a (E F ) 

E F — Eio 



log 



+ g(x) 



Now we define a reference value for E^ 10 , u ^io", by the 
demand that the value of the second line is zero for Eio = 
Eiq. The integral (13) has then the form 

ttiRTA ( ^10 — ^10 k B T 



aE F ' a(E F ) J ' 
where the definition of F RTA is 



{L(\E F \ + \x\k B T )+L>(E F -\x\k B T )}. F RTA (a,b) = J dx + £ 



)) 2 + K, b (x) 



This is integrable with respect to y over the interval be- 
tween one and infinity and independent of the tempera- 
ture between zero and Tq. With Lebesgue we can con- 
clude that the integral of the function (16) of y really goes 
to zero. The terms with a minus in front of y we can treat 
in the same way. For this we introduce functions L < (e) 
and T < (s) in analogy to the above method and demand 
the corresponding property of b(e) of assumption (17). 



with 



t> a ,b( x ) = 9{ x ) ~ CL + xb + log(b). 



(18) 



The corresponding integral in the formula for the linear 
conductance within the DSO in the infinite-U case can be 
represented in an analogous way. The difference is that 
in the numerator we get — f'(x)(l + /(#)) instead of 
only -7r 2 f(x). 
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